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On) ' Abstract 

' I review X-ray observations of accretion-powered millisecond pulsars and current 

theories for formation of their spectra and pulse profiles. 
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O ' 1 Introduction 



in 
o 

^ I A number of rapidly spinning neutron stars in low-mass X-ray binaries were 

Qh' discovered with the Rossi X-ray Timing Explorer [RXTE) in the recent years. 

Q ■ These discoveries confirm the ideas on the formation of radio (recycled) mi l- 

lisecond pulsars in low-mass X-ray binaries (see review by lBhattacharvalll995[ ) . 
! Accretion of matter onto a neutron star results in its spin increase to millisec- 

ond periods if the magnetic field of the star is below about 10^ G. 



Thirteen sources show nearly coherent oscillations for a few seconds during X 



ray b ursts at frequencies ranging from 270 to 619 Hz (see lStrohmaver &: Bildstenl . 



2nn.i for a review). These are now called nuclear-powered millisecond pul- 



sars. The number of accretion-powered millisecond pulsar (AMSP) showing 
pulsations in the persistent emission reached seven by June 20 5. A m ore ob- 
servationally inclined review of AMSP is given bv IWiinandsl Here I 



concentrate on the results of the X-ray spectroscopy, analysis of pulse profiles, 
and our present theoretical understanding. 



2 Pulsars parameters 



The first real AMSP SAX J1808.4-3658 was discovered in 1998 by (|Wiinands fc van der Kli^ . 



1998f ). Now (September 2005) there are seven AMSPs with spin frequencies 
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Table 1 

Parameters of the accretion-powered millisecond pulsars 





Source 


p a 

orb 


u ■ ^ 

spin 


ax sin 


f ^ 

J X 




M ■ ^ 






(min) 


(Hz) 


(It-ms) 


(Mq) 




(Mq) 


1 


SAX J1808.4-3658 


121 


401 


62.809 


3.779 X 10 


-5 


0.043 


2 


XTE J1751-305 


42.4 


435 


10.113 


1.278 X 10 


-6 


0.014 


3 


XTE J0929-314 


43.6 


185 


6.290 


2.9 X 10" 


7 


0.0083 


4 


XTE J1807-294 


40.1 


191 


4.75 


1.49 X 10" 


7 


0.0066 


5 


XTE J1814-338 


257 


314 


390.3 


2.016 X 10 


-3 


0.17 


6 


IGR J00291+5934 


147 


599 


64.993 


2.813 X 10 


-5 


0.039 


7 


HETE J1900. 1-2455 


83.3 


377 


18.39 


2.00 X 10" 


6 


0.016 



" Orbital period; ^ neutron star spin frequency; projected semimajor axis; pulsar 
mass function; ^ minimum companion ma s s for a My = 1.4Mf7^ neutron sta r. 
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from 185 Hz up to 599 Hz (see Table 1). The fastest AMSP, IGR J0029 1+5934 
with the period of just 1.67 ms is the fifth fastest among all known pulsars 
(including radio- and nuclear-powered MSPs). The last AMSP HETE J1900.1- 
2455 was discovered in June 2005. 



AMSPs show pulse frequency variations. These observations are very impor- 
tant for understanding of the evolution of the neutron stars in low-mass X-ray 
binaries towards radio MSPs. They would also shed some light on a compli- 
cated problem of the interaction of the magnetosphere with the accretion flow. 
One expects a spin-up rate 



V 



3.7 X 10 



-13 



^37 



45 



1/2 



M 



2/3 



1.4M0 



spin 

600" 



-1/3 



Hz s~ 



where (notation Q = lO^Q^^ in cgs units is used) I is the neutron star mo- 
ment of inertia, L is the luminosity, rj is the accretion efficiency, and i?co 
are the magnetospheric and corotation radii. Some reported ;> are, however 



negative implying pulsar s lo wing down d uring the outburst (|Gallowav et al 



2OO2I : iMorgan et all . l2003l ). iMarkwardtl (|2004l ) finds wild swings in the ap- 
parent spin frequency of both signs resulting in the total fractional phase 
shift less than 0.15. These could be, however, artifacts of the pulse pro- 
file variations (see Sect. 4). The reported positive z> ~ 8 x 10"^'^ Hz/s for 
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IGR J00291+5934 ()Falanga et al. . 2005bl ) is larger than expected by a factor 
of 5 (since L37 ~ 0.37). A positional error of ~ O'.'T could result i n such a large 



deriv ative, while the radio position is known with a C'l accuracy ()Rupen et al. 



2004 ). The motion of the radio source itself does not produce an error larger 
than ~ 0'.'3 (for the distance of 5 kpc) a week after the outburst, confirming 
the reality of the pulsar spin-up and possibly implying a small moment of 
inertia /. 



Accretion-powered pulsars reside in very compact binary systems with orbital 
periods ranging from 40 min to 4.3 h. Surprisingly three out of 7 pulsars 
have 42 ± 2 min orbits. Adding to this set also 4U1626-67, 4U1916-05 and 
X1832-330 (in globular cluster NGC 6652) with orbital periods of 42, 50, 
and 44 min, respecti vely, it becomes clear that t his interesting fact deserves 
some explanation fsee lNelson fc R,aPDaDortLl2003L for a possible scenario). The 
pulsar mass function 



= (M,sinz)V(M, + = An^a.simf/GP, 



j2 
orb 



is very low for all the objects implying extremely low companion masses con- 
sistent with degenerate white (helium or carbon-oxygen) or brown dwarfs 
(iBildste n fc Cha krabartv'.'200l':'Markwardt et al1 . l2002l : lGallowav et al.Ll2002l : 
IPalanga et all l2005ai .Gallowav et al.. . .200,^1 except XTE J1814-338 which 
contains probably a hydrogen- rich star ( Krauss et al. . 2005[ l. All AMSPs are 
transients with the outbursts repeating every few years and lasting a few 
weeks. They have a rather low time-average accretion rate of ~ 10~^^Mq/jt, 
which could be the main reas on for magnetic field s till to be strong enough 



for pulsations to be observed ( Gumming et al. . 200]| ). 



3 Broad-band X-ray spectra 



The broad-band coverage of the RXTE together with XMM gave a possibility 
to study the spectra of AMSP in great details. The spectra can be modelled 
by three components: two soft, thermally looking ones below a few keV and a 
power-law like tail (see Fig. 1). 

The two soft components which can be modelled as thermal emission from 
a colder {kT ~ 0.4 — 0.6 keV) accretion disc and a hotter (~ 1 keV) spot 
on the neutron star surface. The softer comp onents in XTE J1751-305 an d 
XTE J1807-2 9 4 are s tudied with the XMM bv lGierhhski fc PoutanenI (|2005); 



Falanga et al.l ( 2005a ). The obtained inner disk radius i?in ~ 10—15 km/ Vcos? 
is consistent with the flow disrupted by the neutron star magnetosphere within 
a couple of stellar radii. The hotter black body normalization corresponding 



to the area of ~ 30 — 100 km^ and its pulsation are consistent with it being 
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Fig. 1. Spectrum of XTE J1751-305 as observed by XMM/EPIC-pn, RXTE/PCA 
and HEXTE. The models consists of the multicolour disc (peaking at ~ 1 keV), 
single-temperature blackbody (peaking at 3 — 4 keV), and thermal Comptonization 
of the blackbody photons. The dotted curves show unabsorbed spectral components. 
Panel (a) shows a model where seed photons for Comptonization have the same 
temperature as the blackbody photons, while ii i panel (b) the seed photons are 
hotter. The lower panels shows the residuals. See Gierliiiski fc PoutanenI (|2flnfj ) for 
details. 

produces in a spot at the neutron star surface. 



A power-law tail (having spectral photon index F ~ 1.8 — 2.1) shows a cutoff 
around 100 keV and can be fitted with thermal Comptonization. The electron 
temperature of the Comptonizing medium is around kTe ~20-60 keV and 
Thomson optical depth of tt ~0.7-2.5 (for a plane - parallel slab geometry) 
f Gierlinsj<:i et a L, 20021: Poutanen fc Gierlinski . 120031 : iGierlinski fc Poutanen . 



200,4 iFalanga et al 



2005a|). Weakness of the Compton reflection from the 



disk indicates that the solid angle covered by the disk as viewed from the main 
emission source (accretion sho ck) is small, being consistent with R^^ ~ 40 km 
fjOierUhski fc PoutanenI . \2QQ^ . 



The broad-band X-ray spectra of AMSPs are very similar to e ach other. They 



also show very little variability during the outbursts (see e.g. iGilfanov et al 



19981: IGierlinski fc PoutanenI . l200,4 iFalanga et alibOOSbl ). When fitting spec- 



tra with thermal Comptonization models, one also finds that the product of 
the electron temperature and optical depth is almost invariant (e.g. [kT^, tt] = 
[60 keV, 0.88] in SAX J1808.4-3658; [33 keV, 1.7] in XTE J 1751-305, [37 keV, 1.7] 
in XTE J1807-294, [49 keV , 1.12] in ICR J00291+5934; see lGierUnski fc Poutanen 
l200,4 IFalanga et al.ll2005J bh. The constancy of the spectral slope can be used 
as an argument that the emission region geometry does not vary much with the 
accretion rate. If the energy dissipation takes place in a hot shock, while the 
cooling of the electrons is determined by the reprocessing of the hard X-ray ra- 
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Fig. 2. Normalized pulse profiles of SAX J1808. 4-3658 at different energies during 
the April 1998 outburst. 



diation at the neutro n star surface [two-phase model, e .g. iHaardt fc Maraschi 
I1993I : Poutanen fc S vensson 1996; Malzac et aL 2001 ). the spectral slope is 
determined by the energy balance in the hot phase and, therefore, by the 
geometry. 



4 Phase-resolved spectroscopy, pulse profiles, and time lags 



The pulse profiles from AMSPs are rather close to sinusoidal with peak-to- 
peak oscillation amplitude A = (Fmax — -Fmin)/(-Fmax + -^min) between 4 (in 
XTE J1751-305) and 12 per cent (in XTE J1814-338). At energies above 50 
keV the INTEGRAL data on IGR J00291+5934 show th e increase of the pulse 
amplitude up to ~ 25 per cent (jPalanga et al. . 2005b[ ). Deviations from the 
sine wave are stronger at higher energies (see Pig. 2). The harmonic content 
also is stronger when A is larger (e.g. the harmoii i c-to-f undamental ratio is 
~ 0.03 in XTE ,117 51-305 (lOierhnski fc Poutanenl . liicih . while it is ~ 0.33 
in XTE J18 14-338 (IStrohmaver et all hOOA 



Pulse profiles at higher energies reach their peaks at an earlier phase rela- 
tive to the soft photons resulting in th e soft time lags. I n SAX J1808. 4-3658 



tive to tne sort pnotons resulting m tn e soit time lags, m aAA jiaug.4-dbos 
(ICui et al.Lll998[lGierlinski eta l.'. '2002 1), XTE J175 1-305 (|GierUnski fc Poutanen 
200.51 ). and XTE .10929-314 (|Gallowav enil . l2002l ) the lags increase (in abso- 
lute value) with energy up to about 7-10 keV after the y saturate (see Fig. 3, 
left panel). In IGR J 0029 1+5934 (see Pig. 3, right panel: iGallowav et al.ll20oi 
Palanga et al 1l200.5b[ ). the behaviour is more complicated and the lags decrease 
between 7 and 15 keV, seemingly saturating at higher energies. 



One can point out that the contribution of the black body decreases exponen- 
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Fig. 3. Time lags a s a function of energy for SAX J 1808.4-3658 (left pan el; 
Gierlinski et~aDl2nn2h and IGR J00291+5934 (right panel: iFalanga et alJEnoHbh . 
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Fig. 4. Pulsation of the two model com ponents, blackbody and Comptonization, in 
XTE J1751-305. From lOierlinski fc P outanen (2005). 



tially with energy and the lags increase significantly at the same time. Above 
~ 7 keV, the blackbody's contribution is negligible and the lags saturate. 
When fitting phase-resolved spectra with a two-component model (blackbody 
+ Comptonization), one clearly sees that the normalizations of the compo- 
nents do not vary in phase, with black body lagging the Comptonized e mis- 
sion (see Fig. 4 and iGierlihski etHI 120021 : iGierlihski fc PoutanenI 1200,51 ). In 
SAX J1808. 4-3658 the profile corresponding to the Comptonized emission is 
clearly non-sinusoidal. Difference in profiles of the two components can b e ex- 
plained by their different emission pattern Poutanen fc Gierlinski ( 2003[ l. Its 
variation may cause pulse shape change resulting in the swings of i> discussed 
in Sect. 2. 



5 Modelling the pulse profiles of AMSPs 



One believes that the bulk of the X-ray emission observed from AMSPs orig- 
inates from polar caps where the gas stream channels by the neutron star 
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Fig. 5. Pulse profile for a small black body spot at the neutron star surface, (a) Dot- 
ted curve shows the profile for a slowly rotating star in Newtonian approximation. 
Gravitational light bending reduces variability amplitude (dashed curve). Doppler 
effect due to rapid rotation skews the profile (dot-dashed curves). Accounting for 
the light travel delays further modifies the profile slightly (solid curve), (b) Doppler 
factor as a function of phase. We assumed neutron star mass M = 1.4M0 and radius 
R = 10.3 km, rotational frequency u = 600 Hz, the inclination i = 45°, and the 
polar angle of the spot center 6 = 45°. 

magnetic field impacts the stellar surface forming a shock. This is supported 
by a weak energy dependence of the variability amplitude and a fairly con- 
stant spectral shape as a function of pulse phase (any additional source of 
radiation would have to have a spectrum identical to that of the shock), as 
well as the observed broadening of the pulse peak in the power-spectrum of 
SAX J1 808.4-3658 due to m odulation of the aperiodic variability by the spin 



period ()Menna et all 120031 ). 



5.1 Computing pulse profiles 



Pulse profile shape and variability amplitude carry the information about 
the compactness of the neutron star, spot size and its position at the star, 
and the emissio n pattern. Gravitational light bending tends to decrease the 
pulse amplitude (IPec henic k et al Lll983l: IPiffert fc Meszarosl . Il988l : IPagel . Il99,4 
Leahv fc Li . 199^ BeloborodoX 2002f ) . An e arlier attempt to i ncorp orate the 



effect of rapid rotation using Kerr metric ( Chen &: Shaharnl 1989i) did not 
include the dominant Doppler boost caused by the star's surface mo tion in 
the locally non- rotating frame. This was criticized by Braie et al. ( 2000l ). They 
also noticed that Doppler boost and light travel time delays are important for 
shaping the profiles, while the frame dragging affe c ts the m at a ~ 1% level. 
The Monte-Carlo simulations used by iBraie et al. I (|2000h for the light curve 
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calculation are, however, extremely time-consuming and are hardly suitable 
for data fitting. 



As frame dragging is not so important, one can consider Schwarzschild metric 
instead. One can account for Doppler effect making Lorentz transformation 
from the frame rotating with the star to a non-rotating frame and then follow 
light trajectory in the Schwarzschild space-time to infinity ("Schwarzschild + 
Doppler" or SD appro ach) . An attempt to formulate the problem was done by 
Miller fc LambI ( IQQSi l who considered the spot and the observer in the equa- 



torial plane only. Oscillation waveforr ns and amplitud e s for a rbitr ary spot and 
observ er position were computed by IWeinberg et al.l (|2001l l and iMuno et al. 
(120021 ). These paper provided no details on how to account for effect of aberra- 
tion on intensity and projected spot area. A formalism for computing the pulse 
profiles accounting for Doppler boosting, rela tivistic aberration and gravita- 
tional bending appeared for the first time in Poutanen &: Gierlihski (j2003|). 
Effects of the anisotropy of Comptonized radiation from a shoc k on th e light 
curves and polarization were studied by Viironen fc PoutanenI (12004 ) . Ana - 
lytical formulae for oscillation amplitudes are presented bv lPoutanenI ( 2004 ). 
Light curves from realistic sp ots produced by accret i on on to inclined mag- 
ne tic dipole are computed by Kulk arni fc Romanova^ ('2005') usi ng formalism 
of Poutanen fc Gierlihski ( 2003[ ): IViironen fc Pouta nen (200l: iBeloborodovl 
( 2002| l. Accuracy of the SD approach is discussed by Cadeau et al. ( 2005( l. 



5.2 Main effects 



Let us now briefly describe the main effects shaping the pulse profiles of AM- 
SPs. Without general or special relativistic effects, a small black body spot 
would produce sinusoidal variations (with possible eclipses) due to a change of 
the projected area. Light bending reduces the variability amplitude (compare 
dotted and dashed c urves in Fig. 5), while the pulse remains almost sinusoidal 
()BeloborodovLl2002h . Relativistic aberration and Doppler boosting modify the 
observed flux for rapidly spinning star. The projected area is changed by the 
Doppler factor 5 due to aberr ation, and the frequency -integrated specific in- 
tensity is multiplied by 5 (see iPoutanen fc Gierlihskil .12003). Since 5 reaches 



the maximum a quarter of the period earlier than the projected area (Fig. 5b), 
the pulse becomes skewed to the left (see Fig. 5 and compare it to Fig. 2). 
Light travel time delays slightly modify the profile further (note, that time 
delays are already accounted for in the flux calculations by one of the Doppler 
factors) . 
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5.3 Oscillation amplitudes 



Due to the combined action of the Doppler effect and time delays, the pulse 
profile deviates from a simple sine wave. (We consider below the situation 
where the spot is always visible.) The relative amplitudes (in the bolometric 
signal) of the harmonic Ai to that of th e fundamental A n grows approximately 
linearly with the rotational frequency (|Poutanenl . 12004^ : 



An 



5 2'kR 

2 ~ir 



fspm sin i sin 6 



0.16 



R 



spm 



10 km 300 Hz 



sini sin 6', 



where i is the inclination and 9 is the polar angle of the spot center. Even a 
slight deviation from the isotropic emission (e.g. in the form I {a) oc 1 + 6 cos a, 
where a is the angle from the spot normal and h is the anisotropy parameter) 
leads to a high harmonic conteri t (|Poutanen fc Gierhhskil . 120031 : iPoutanenl . 
2004 IViironen fc Poutanenl . l2004[) : 



Ai _ h{l-RjR)/2 sin i sin ^ 

~ l + 2h[RjR+{l- RJR) cos i cos 6] ' 



where i?s = 2GM/c^. The h armonic content A^ /An as well as the total r r as are 
proportional to sin z sin e fiBeloborodovL 120021 : IPoutanen fc Gierhhskil 1200.1 
2004 IPoutanenl . 120041 . Therefore, one expects harmonics to be relatively 
stronger when rms is larger which is exactly what is observed. Increasing 
the spot size leads to a reduction of the amplitude of the ha rmonic first and 
of the total rms later (jWeinberg et al.l . l200ll : IPoutanenl . 120041) . 



If the spectrum has a sharp cutoff, the rms amplitude of the pulse at en- 
ergies above the cutoff incre a ses dramatically. The flux v aries as oc 
f Poutanen fc Gierlihski . 20031 : Viironen fc Poutanen . 2004), where T(E) = 
1 — din Fe/ din E is the energy-dependent photon spectral index, resulting 
in the rms energy dependence oc 3 + T{E). For the spectrum in the form 
Fe oc E~^^°~^^ exp{—[E / Ec]'^), which is typical for Comptonization, the pho- 
ton index T{E) = Tq + ^^E/E^f. At low energies, F ^ Fq, and rms is a 
very weak function of energy. Close to the cutoff, the spectral index rapidly 
increases resulting in a significant grow of the rm s , as ob served in the INTE- 
GRAL data for IGR ,100291+5934 (IFalanga et al.Ll2005bl ). In case of coherent 
oscillations observed during X-ray bursts, a linear increase of the rms with 
energy is expected, since t he spectrum is close to a blackbody wh ich has an 
exponential cutoff (see also Miller fc LambI 19981: Muno et al. 20031 ) . 
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Fig. 6. Constraints on the neutron star mass-radius relation obtained by fitting the 
pulse profiles of SAX J1808. 4-3658 (filled circles with erro r bars) together with a 
set of equations of state for neutron and strange stars. From Poutanen Sz Gierliiiskil 



5.4 Constraints on the neutron star equation of state and spot parameters 



The amphtude of variability and the pulse shape (harmonic content) can be 
used to put constraints on the compactness of the neutron star, position, size 
and the emission pattern of a radiating spot. In case of the coherent oscilla- 
tions observed during the X-ray bursts the constraints on the compactness are 
not v ery tight (jMiller bambl . Il998l : iNath et abl 120021 : iRhattacharvva et a1 " 



2005| ). A very weak harmonic content led to the conclusion that the bright re- 
gion (or the line of sight) must be near the rotational pole or must cover half 
of the star (IMuuo et al.[ l2002[ ). Fitting the pulse profile shape (but not the 
amplitude) of the osci llations observed during the X-ray burst of XTE J1814- 
338 was performed bv tehattacharyva et al. ( 2005h . who obtained the preferred 
spot polar angle 6 = 90° ±30°. This is highly improbable si nce an identical (in 



shape and phase) pulse observed in the persistent emission (jStrohmaver et al 



2OO3I ) would require the magnetic poles to lie at the equator and the antipodal 
spot should then be visible too. 



For AMSPs the statistics in the pulse profiles is much better since the pulse is 
folded over a longer observational period (days or weeks rather than seconds 
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as in the case of X-ray bursts). This allowed to obtain tight constraints on the 
neutron star compactness in SAX J1808. 4-3658 as well as the emission pattern 
(consistent with Comptoniza tion from the plane- parallel shock), inclination, 



and the spot polar angle fsee iPoutanen fc GierlihskL.2003» . and Fig. 6). 



6 Summary 



We have reviewed here X-ray observations of the AMSPs concentrating on the 
broad-band spectra and the energy-dependent pulse profiles. Recent advances 
in the theoretical modelling of the spectral and timing characteristics are also 
described. At the observational side, we are missing the details of the pulse 
profile and spectral evolution during the outbursts. One would expect that 
the changing accretion rate would force a change in the accretion pattern on 
to a neutron star surface, causing corresponding profile variations. Detailed 
studies of the reflection amplitude and frequencies of quasi-periodic oscillations 
should reveal changes in the inner disk radius during outbursts. The pulse 
profile should also react correspondingly as the visibility of the antipodal spot 
depends on the disk radius. Combined with the spin-up rate observations, 
these would shed some light on a complicated problem of the magnetosphere 
interaction with the accretion flow. On the theoretical side, we need a detailed 
model describing the dynamics of the accretion flow onto the neutron star 
surface including radiation feedback. It is not clear yet what fraction of the 
material follows the magnetic field lines towards the poles. The impact of the 
realistic spot geometry and the radiation pattern on the pulse profiles and the 
derived parameters should be studied further. 
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